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ABSTRACT 

We apply a wind model, driven by combined cosmic-ray and thermal-gas pressure, to the Milky Way, and 
show that the observed Galactic diffuse soft X-ray emission can be better explained by a wind than by previous 
static gas models. We find that cosmic-ray pressure is essential to driving the observed wind. Having thus 
defined a "best-fit" model for a Galactic wind, we explore variations in the base parameters and show how 
the wind's properties vary with changes in gas pressure, cosmic-ray pressure and density. We demonstrate the 
importance of cosmic rays in launching winds, and the effect cosmic rays have on wind dynamics. In addition, 
this model adds support to the hypothesis of Breitschwerdt and collaborators that such a wind may help explain 
the relatively small gradient observed in 7-ray emission as a function of galactocentric radius. 
Subject headings: ISM:outflows - ISMxosmic rays - ISM:magnetic fields - Galaxy : evolution - X-rays:diffuse 
background 



1. INTRODUCTION 

Large-scale galactic outflows are usually considered in 
the context of starburst galaxies or Active Galactic Nuclei 
(Veilleux et al. 2005). These outflows are interesting not only 
intrinsically (what drives the outflow?) but for the interstellar 
and intergalactic media (how is the host galaxy affected, and 
what metals are ejected from the galaxy?). 

To examine these questions, we have built a thermal and 
cosmic -ray driven wind model. Our investigation into such 
models was first inspired by observational hints that the Milky 
Way may possess a kiloparsec-scale wind; this paper further 
explores that possibility. To motivate this study, we first intro- 
duce the observational evidence for a Galactic wind (§§L1 
and L2) and then introduce the cosmic -ray and thermally 
driven wind model (§2). In §3, we calculate the X-ray emis- 
sion from this wind, and then compare it to the observations. 
After finding the best-fit wind model, we then explore the 
parameter space around that model (§4) to understand more 
about how the wind is modified by varying the input and fit 
parameters. Our conclusions are given in §5. 

But first, observational hints for an outflow from our own 
Galaxy. 

1.1. X-ray Observations 

The Milky Way does exhibit clues that it might drive a 
large-scale wind. The first of these is an enhancement in 
the diffuse soft X-ray emission, stretching over the longi- 
tude range -20° < / < 35° with an emission scale height in 
the southern Galactic hemisphere of b ^ -17° (see Fig. 1). 
This emission was first noted by Snowden et al. (1995), who 
modeled it with an isothermal plasma with a temperature T = 
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4 X 10^ K, a midplane electron density of ne.midpiane 3.5 x 
lO"-' cm"-', and a midplane thermal pressure of Pg.midpiane/^ ^ 
2.8 X 10"^ cm"-' K. At approximately the same time, Bre- 
itschwerdt & Schmutzler (1994) suggested that the average 
all-sky X-ray emission (not only that emission in the region 
defined above) in all ROSAT bands might be explained by 
delayed-recombination in a large-scale cosmic -ray and ther- 
mally driven wind (see also Breitschwerdt & Schmutzler 
1999). 

Later, Almy et al. (2000) used intervening absorption to 
show that at least half of the central, enhanced X-ray emis- 
sion lies more than 2 kpc from the sun (see also Park et al. 
1997, 1998). Since that measurement was made in the Galac- 
tic plane, where the absorption is strongest, it was inferred 
that most of the emission observed at higher latitudes lies be- 
yond that 2 kpc distance. Almy et al. (2000) also improved 
on previous modeling efforts: that work presents a model of 
the emission due to a static poly tropic gas (with 7 = 5/3), 
and very importantly, includes the effects of known back- 
ground components, such as the stellar background, extra- 
galactic background, and an additional isotropic background 
(to fit high-latitude emission). For comparison, their model 
had a central temperature of To = 8.2 x lO*" K, a central elec- 
tron density of We = 1 • 1 x 10"^ cm"^, and a central pressure of 
Pg,o=1.8x 10^ cm"3 K. 

We will compare our results with this static polytrope 
model to investigate whether a wind model for this emission 
is feasible. 

1 .2. Cosmic Ray Source Density 

Another indicator of a Galactic wind comes from measure- 
ments of the density of cosmic rays as a function of Galacto- 
centric radius, R. The source density of cosmic rays can be 
determined via 7-ray emission: the production of 7-ray pho- 
tons with energies exceeding about 50 MeV is dominated by 
collisions of cosmic rays with gas in the interstellar medium 
(Bloemen et al. 1984). Since the galaxy is largely transpar- 
ent to such high-energy photons, the 7-ray emissivity at those 
energies yields the cosmic ray source density. 

If cosmic rays are produced in supernovae remnants, then 
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Fig. 1 . — X-ray emission at 3/4 keV (tlie "R45 band") as seen by ROSAT (Snowden et al. 1997). These observations suggest a "Galactic X-ray Bulge", seen 
most clearly in the southern Galactic Hemisphere, and stretching over the Galactic longitude range, /, from |/| < 30° and down to approximately -15° in Galactic 
latitude. This paper asks whether the X-ray bulge in the southern Galactic Hemisphere can be explained with a combined thermal and cosmic-ray driven wind. 



since the source density of supernovae remnants seems to 
increase with decreasing R, the cosmic -ray source density 
should increase as well. However, it has been known for some 
time (e.g, Bloemen 1989) that the inferred cosmic -ray source 
density is relatively flat, compared to the supernova density, 
as a function of R. 

There has, however, been some debate about whether su- 
pernovae remnants are an accurate tracer (since those surveys 
are subject to various selection effects; see, e.g.. Strong et al. 
2004, and references therein). Recent surveys of the pulsar 
population (Lorimer et al. 2006) also show that the pulsar 
source density increases towards the center of the Galaxy, as 
shown in Figure 2. This is true irrespective of the model of 
how «e varies in the disk, although the magnitude of the pulsar 
population gradient with R depends strongly on the model. 
So, there remains a mismatch between the observed source 
density of cosmic ray "producers" and the cosmic rays them- 
selves. 

It has already been pointed out that the observed slow rise 
in cosmic rays may be due to a wind emerging from the disk, 
advecting cosmic rays outwards (Bloemen et al. 1993; Bre- 
itschwerdt et al. 2002). In the case of Bloemen et al. (1993), a 
wind model was applied to the entire Galactic disk; as a result, 
only a very slow wind was found to be compatible with the in- 
ferred cosmic -ray source density. In contrast, Breitschwerdt 
et al. (2002) applied their cosmic-ray and thermally driven 
wind model, where the wind velocity varied as a function of 
radius and height; they also took into account anisotropic dif- 
fusion. With this model, a small radial gradient in the cosmic 
ray source density could be explained. 

An alternate explanation for this slow change in the cosmic 
ray population with R was proposed by Strong et al. (2004), 
who found that a radial variation in the Wco-lo-N{H2) ratio by 
a factor of 5 to 10 could explain the 7-ray observations. In this 
paper we primarily address the question of the origin of the 
diffuse, soft X-ray background emission; we will, however, 
concentrate on a large-scale wind model, keeping in mind its 
possible application to the cosmic ray source density. 

2. A COSMIC RAY AND THERMALLY DRIVEN WIND 
MODEL 
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Fig . 2. — Comparison of two different calculations of the pulsar population 
as a function of Galactocentric radius (Lorimer et al. 2006) vs. the cosmic 
ray source density implied from the observed 7-ray emissivity (Strong et al. 
2004). The two different curves for the pulsar distribution result from as- 
suming a smooth distribution of in the Galaxy (Lyne et al. 1985), or a 
clumped distribution, using Cordes & Lazio (2002) and Faucher-Giguere & 
Kaspi (2006); for details, see Lorimer et al. (2006). The fact that the cosmic- 
ray distribution does not seem to follow the pulsar population has been known 
for some time (Bloemen 1989), but there is no consensus on the reason. A 
cosmic-ray and thermal pressure-driven wind may help explain the cosmic- 
ray source population. 



To understand the observations outlined above, we must ad- 
dress both the thermal (soft X-ray) gas and the cosmic ray 
source density. Therefore, in investigating the possibility of 
an outflow, we require a Galactic wind model that includes 
cosmic -ray pressure (exerted via Alfven waves) and thermal 
gas pressure. We outline such a model in this section. 

Our chief motivation here is to fit the wind model to the 
observed large-scale soft X-ray emission, and as such, we re- 
quire a relatively simple model that can be computed quickly 
to compare with the observations, but one that takes into ac- 
count the physics of cosmic-ray interactions with the ther- 
mal gas. In addition, we are interested in building intuition 
into the differences between pure thermaUy-driven winds and 
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winds with a significant cosmic -ray component (see § 4). We 
are not the first to address this; a very suitable model has al- 
ready been developed by Breitschwerdt et al. (1991, hereafter 
BMV91) and further advanced in later papers (Breitschwerdt 
et al. 1993; ZirakashviU et al. 1996; Ptuskin et al. 1997; Bre- 
itschwerdt & Schmutzler 1999; Breitschwerdt et al. 2002). 
This work had built on previous analyses of the possibility of 
cosmic-ray driven winds (Ipavich 1975; Breitschwerdt et al. 
1987). In what follows, we explain this ID, semi-analytic 
wind model in detail, with particular attention to the differ- 
ences between our model and that of BMV91 . 

This ID model must assume a particular geometric cross- 
section for the wind; we adopt a flared-cylinder geometry 
shown in Figure 3, as was adopted in BMV91 and used in 
studies of coronal holes on the Sun (e.g., Kopp & Holzer 
1976). In this geometry, the wind flows along a tube of ap- 
proximately constant cylindrical cross section up to a height 
z ~ Zbreak, after which the area increases as z": 



A(z) =Ao 



1 + 



Zbrtzk 



(1) 




midplane 



Fig. 3. — The geometry of a flow tube for the ID wmd model. The wind 
starts at the Galactic midplane with cross-sectional area Aq, which is roughly 
constant up to z ~ Zbreak, after which A oc z™ . 



We envision that this geometry is set by approximate pres- 
sure confinement in the plane of the galaxy up to Zbreak; above 
that height, the area diverges spherically (thus a = 2). Unlike 
BMV91, we allow Zbreak to change to find the best-fit wind 
model. We do not 'freeze' this parameter, as we have no a pri- 
ori value for Zbi-eak for this particular wind geometry (which, 
in contrast to BMV91, covers only a part of the Galactic disk, 
and so does not require a scale-height comparable to the scale- 
length of the star-formation disk for an approximately spheri- 
cal geometry, for instance). Thus, Zbreak is a fit parameter for 
our model, in the context of the simplified geometry that we 
are assuming. Also unhke BMV91, we assume that the wind 
is launched from the Galactic midplane; this wiU be addressed 
in detail when the interaction of cosmic rays, magnetic fields, 
and thermal gas is discussed, below. 

With this geometry, and assuming no additional mass load- 
ing of the outflow, the wind has a simple equation of mass 
conservation: 

^(pvA) = 0, (2) 
dz 

where z is the height above the Galactic midplane, p is the 
density in the wind and v is the wind velocity. 

Next, how is the wind driven along this flow tube? We 
wish to consider the possibility that cosmic-ray pressure, 
Alfven wave pressure, and thermal pressure are aU impor- 
tant components in driving a Galactic wind. At first glance, 
this may not seem productive, as cosmic rays seldom inter- 
act directly with any particle in the galaxy: the probability of 
any cosmic ray particle colliding with matter in the Galaxy in 



their lifetime is of order 10~^ (see, e.g., Kulsrud 2005). How- 
ever, cosmic rays are observed to have a very small anisotropy 
(about 1 part in lO'*), which seems at odds with this low col- 
lision rate. This small anisotropy is explained by pitch-angle 
scattering of cosmic rays by Alfven waves in the ISM. It has 
been shown that the cosmic rays which supply most of the 
pressure {E < 100 GeV) can generate these waves via the 
"streaming instability" (e.g., Wentzel 1968; Kulsrud & Pearce 
1969; Kulsrud & Cesarsky 1971). This instability amplifies 
waves with wavelength of order the cosmic ray gyroradius 
when the bulk velocity of cosmic rays along the fieldlines ex- 
ceeds the local Alfven speed. If the cosmic ray mean free path 
is much shorter than global lengthscales in the problem, the 
cosmic rays can be described as a fluid which moves down 
its pressure gradient at velocity va relative to the thermal gas, 
while transferring momentum and energy to the waves, which 
in turn transmits them to the thermal gas (e.g., Skilling 1975; 
Drury & Volk 1981; McKenzie & Webb 1984). We adopt fliat 
picture here. 

In computing the interaction of cosmic rays, Alfven waves, 
and thermal gas, we treat the cosmic rays as an ultra- 
relativistic polytropic gas with 7^ = 4/3. (This is an approx- 
imation, as 7c is not exactly 4/3; see Ensslin et al., 2006.) 
Meanwhile, we treat the thermal gas as having a polytropic 
index, 7^, of 5/3. Much as in BMV91, we derive equations 
for the change in gas and cosmic-ray pressure with height in 
the wind: 



dz 



-7f= k -7c(7.-l) 



Pc 1 Ma + \ 
p MaMa + 1 



and 



dPc _-fcPcMA + jdp 



^, (3) 



(4) 



dz p Ma + 1 dz' 

where Pg and Pc represents the gas and cosmic-ray pressure, 
Cg gives the sound speed in the gas, and Ma = v/va is the 
Alfven Mach number for the wind. In the gas pressure equa- 
tion. Equation 3, the first term in parentheses simply relates 
the change in gas pressure to the change in density as the gas 
accelerates and expands in the flow tube. 

The second term in the parenthesis of Equation 3 represents 
the coupling of cosmic -ray generated Alfven waves to the gas; 
that term gives the heat input to the gas from the damping 
of those cosmic-ray generated waves. These waves represent 
the dominant coupling between the cosmic rays and the ther- 
mal gas; this process heats the gas despite the drop in density 
with height (hence the negative sign for this term). As the 
cosmic-ray generated Alfven waves are immediately damped, 
they do not add to the wave pressure, and hence we do not 
follow their evolution. In our models, the wave pressure at 
the base of the wind is set to zero, reflecting the smaU wave 
energy density in the Galactic plane (6B/B ~ 10"^ from Kul- 
srud & Pearce, 1969, where SB represents the Alfven wave 
perturbation to the large-scale magnetic field, B). Any energy 
transfered from the cosmic rays to the waves is immediately 
input to the gas, as in Equation 3 above. The inclusion of this 
immediate wave damping in all of the models presented here 
is an important difference between this work and most of the 
models in BMV91; we also note that wave damping was con- 
sidered in much more detail in the later papers of ZirakashvUi 
et al. (1996) and Ptuskin et al. (1997). 

This immediate damping of the cosmic-ray generated 
Alfven waves is important for two reasons. First, left 
unchecked the wave pressure can easily grow to such mag- 
nitudes that the ratio of the perturbed magnetic field to the 
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large-scale magnetic field, SB/B, exceeds unity. In this limit, 
the derivation of the above equations becomes suspect, as the 
system becomes nonlinear. 

The second objection to unlimited Alfven wave growth is 
the presence of rapid damping mechanisms in the ISM. In 
particular, non-linear Landau damping (Kulsrud 2005) will 
quickly remove energy from the waves. The rate of non-Unear 
Landau damping is of order 



7NL ■ 



SBvA^ 



(5) 



where fi is the ratio of the gas pressure to magnetic pressure, 
(87rPg)/B^, rti is the ion cyclotron frequency, and v, is the 
ion thermal velocity. For typical values for our wind mod- 
els at z ~ 2 kpc, (B = 7.2/iG, 5B/B ~ 10"^), we find tnl ~ 
2 X 10"^ s"', or Tnl ^ 5 x 10^ sec = 1.6 yr. We can compare 
this to the advection timescale, 1 / ((v + v^) • S/Pc/Pc), which 
for a typical wind model is of order lO^'' s, or 4.3 x 10* years. 
From this comparison, we can see that the damping is local, 
since the timescale for non-linear Landau damping is much 
smaller than the advection timescale. Therefore, the local 
damping not only allows the quasihnear equations to be ap- 
pUcable throughout the wind, but is physically quite plausible, 
given the above small damping timescale. 

(We note that these arguments were well-known to BMV9 1 , 
but that the wind models calculated there did not develop large 
6B/B inside the critical point, so such immediate damping 
was not important in those models [Breitschwerdt 2007, per- 
sonal coimnunication].) 

We briefly consider the effect of Alfven waves generated 
by other sources. Cosmic ray streaming is only one source of 
MHD turbulence in the ISM, and, at large scales, by no means 
the dominant one. If the turbulence which couples the cos- 
mic rays to the gas were not primarily due to the streaming 
instability, the model would be substantially modified. The 
Alfven speed which appears in the model should be thought 
of as the mean velocity of the waves in the rest frame of the 
thermal gas; if the waves were isotropic, this speed would be 
zero. The cosmic rays would simply behave as a polytropic 
fluid with adiabatic index jc and there would be no frictional 
heating of the thermal gas. This is not a significant issue, 
however, because of the anisotropy in the turbulent cascade 
seen in both experiments and theoretical studies of MHD tur- 
bulence (see, e.g. Shebalin et al. 1983; Goldreich & Sridhar 
1995; Cho & Vishniac 2000; Milano et al. 2001), which has 
also been invoked in studies of turbulence in supernovae as 
well (Ptuskin & Zirakashvili 2003). The turbulence is ini- 
tially presumably excited at scales of several parsecs or more, 
far above the cosmic ray gyroradius scale of 10'^- 10^^ cm, 
and cascades down to the gyroradius scale through nonUn- 
ear wave-wave interactions. In order to interact, these waves 
must be oppositely directed along the magnetic field. By mo- 
mentum and energy conservation, interactions between such 
oppositely-moving waves yields resultant waves in which the 
component of momentum along the magnetic field Une has 
not increased. However, the perpendicular component can in- 
crease as a result of the interaction, yielding an anisotropy in 
k-space (Shebalin et al. 1983). Thus, at wavenumbers much 
greater than the driving scale, the perpendicular wavenumber 
kj_ much exceeds (Goldreich & Sridhar 1995; Cho & Vish- 
niac 2000). More intuitively, for motions on smaller scales 
and commensurately smaller energies, the turbulent motions 
cannot bend magnetic field lines, and the energy is tranferred 



to motions parallel to the magnetic field, resulting in elon- 
gated eddies (Shebalin et al. 1983; Lazarian 2006). 

So how does such an anistropic eddy affect cosmic ray scat- 
tering? When kj_/k\\ » 1, cosmic rays scatter much less ef- 
ficiently; this is because, when k± » individual cos- 
mic rays sample many uncorrelated perturbations in k± in 
each gyro-period that effectively cancel out (Chandran 2000; 
Lazarian 2006). Thus, the background Galactic turbulent cas- 
cade has little effect on the cosmic rays compared to the waves 
they generate themselves. We note, however, that fast-mode 
waves in the ISM may perhaps be important in scattering cos- 
mic rays (Yan & Lazarian 2003, 2005). For this work, how- 
ever, we retain the model of quickly damped, cosmic-ray gen- 
erated Alfven waves in Equation 3. 

Given the importance of gas heating, one may wonder if 
perhaps a loss term, such as radiative cooling, is also impor- 
tant. We have determined that the radiative cooling is ^ 1 % 
of the total power in the wind, and so is unimportant for the 
wind models presented here. However, we note that for very 
low velocity winds, the dynamical timescale may exceed the 
radiative cooling time, resulting in significant cooling. 

We also ignore cosmic ray diffusion and thermal conduc- 
tivity. Diffusion is important towards the base of the wind 
(Breitschwerdt et al. 1993); thermal conductivity may be im- 
portant there as well: our calculations show conductivity to 
be important for z < 350 pc in the wind (below that height, 
the energy input from conduction dominates adiabatic cooling 
and heating via wave-damping). This is certainly significant, 
and both effects need to be considered in a more detailed wind 
model (see §5). 

Given Equations 2 to 4, and our assumptions about the cou- 
pling of cosmic rays and thermal gas, the above pressure re- 
lations are then coupled together in the wind equation, which 
in its simplest form is 



dv 2 
P'd-z^'*Tz=-P'' 



(6) 



where c*, the "composite sound speed" (see BMV91), is 
given by 

2 _ d(Pg+Pc) 



dp 



(7) 



and where the gravitational acceleration, g, is defined by 
a three-component (bulge, disk, and halo) model given in 
BMV91. We have compared this gravitational potential 
model to the more recent work of Dehnen & Binney (1998), 
and found that, for z > 200 pc, this newer model gives a lower 
gravitational acceleration (by ^ 10% to at most ~ 60%) than 
the simple model of BMV91; closer to the disk, the model 
of Dehnen & Binney (1998) yields a higher gravitational ac- 
celeration by ^ 40%. With the reduced calculational require- 
ments of the simpler potential (especially important near the 
critical point of the wind), and the knowledge that the simpler 
model largely over-estimates the gravitation potential where 
the wind is accelerating, we choose the simpler model as a 
conservative approximation; using the more realistic poten- 
tial, launching a wind from the Milky Way should be some- 
what easier than shown here. 

We solve the wind equation (Eq. 6) in the ID flow tube de- 
fined by the area law (Eq. 1). While integrating, we use the 
simple magnetic flux conservation law, d(BA)/dz= 0. The 
integration is carried out using the dlsode routine in ODE- 
PACK (Hindmarsh 1983). 
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TABLE 1 

Initial (Midplane) Wind Parameters at 
R = 3.5 kpc 



Parameter 


Value 


Fixed? 


Rq Range [Galactocentric] 


1.5 to 4.5 kpc 


Fixed 




1 pc 


Fixed 


no 


1.8 X 10-2 cm-3 


Varied 


Pg,0/kB 


2.6xlff*Kcm-3 


Varied 




2.2x 10* Kcm-3 


Fixed 


Bo 


7.8 /iG 


Fixed 


Zbreak 


4.5 kpc 


Varied 


a 


2.0 


Fixed 



2. 1 . Initial Conditions 

The integrations require that Pg, p, R, and B all be spec- 
ified at the base of the wind. The value of v is not given at 
the base of the wind, as its value is set by requiring that the 
integration pass through the critical point (see § 2.2). Fiducial 
values for these parameters are given in Table 1, and discussed 
in more detail below. 

The Galactocentric launching radius is set by the observed 
geometry of the diffuse X-ray emission. If the X-ray emission 
is roughly centered on Sgr A* (Almy et al. 2000), the extent of 
the emission out to longitudes of ^ 30° implies a scale for the 
outer edge of the wind at /? ~ 4.5 kpc (assuming that the Sun 
is located at ^ 8 kpc). The inner edge of the wind is set by 
the difficulty in launching from the central Galactic potential; 
for ISM-like launching conditions, we have found that models 
with R<1.5 kpc caimot escape the Galaxy. Thus, we set the 
geometry of our wind to be a "thick curtain", launched over 
Galactocentric radii of 1 .5 to 4.5 kpc. To run a fiducial model, 
we set the wind flowtobe at 3.5 kpc. More complicated simu- 
lations (where separate winds are launched over the complete 
range or radii) have been run, and show differences in X-ray 
emission of only ~ 10%. Therefore, the wind models here 
are simply run at ^ 3.5 kpc, and that wind is then applied 
(or "replicated") to the range of wind radii (1.5 to 4.5 kpc) and 
over 27r in azimuth, to fit the observed diffuse X-ray emission. 
This geometry is shown in Figure 4, below. 

This is, of course, a simplification, but one that allows quick 
calculations of the wind's X-ray emission and comparison to 
observations (§ 3), and allows surveys of large areas of param- 
eter space for building physical intuition (§ 4) about galactic 
winds. It is important to note that, of course, at high lati- 
tudes (where z > Zbreak), each tube flares outward; with such 
tubes placed next to each other, they will strongly overlap for 
z ^ Zbreak- This could be a potential problem, but as we will 
see later (in §3.4), the X-ray emission is explained with Zbreak 
being more than a factor of two larger than the emission scale- 
height, so this will not impact our predictions of the X-ray 
emission. It would, however, be important for observational 
tests at higher latitudes. 

We must also choose the initial height from which to launch 
the wind. Again, we differ from BMV91, and choose to 
launch the wind from the midplane of the Galaxy, with zo = 
1 pc. BMV91 choose to launch from zo = 1 kpc due to con- 
cerns about (1) ion-neutral friction due to partially neutral gas 
for z < 1 kpc, and (2) an isotropic Alfven wave field (not gen- 
erated by the streaming instability of the cosmic rays). While 
those are all important considerations for the generic ISM, we 
hypothesize that the hot, ionized medium from which these 
winds are launched is largely free of neutrals, is dominated 
by heating and cosmic rays due to nearby supemovae, and that 



J._l._l J,_.1,_J — l_,J._,L l_.J,_.L midplane 

R,-S R R, R0t 5R 



Fig. 4. — Near the Galactic midplane, the flowtubes are vertical, allowing 
the tubes to be assembled in a simple, large-scale wind model. In this way, 
our ID wind models are assembled into a large-scale wind to compare to the 
observed diffuse soft X-ray emission. 



cosmic -ray generated waves dominate the scattering process 
(see §2). Also, for meaningful comparison with observations, 
we must model the wind below z ~ 1 kpc, as the scale height 
of the diffuse soft X-ray emission is z ~ 2 kpc (Snowden et al. 
1997). 

In order to fit the wind models to the observed diffuse 
Galactic X-ray emission, the variables Pg, p, and Zbi-eak will 
be left to "float" (hence allowing the temperature and emis- 
sivity of the gas the possibihty of matching the observations). 
But for comparison, it is helpful to consider what prior ob- 
servational constraints we can place on these values. Un- 
fortunately, it is difficult to constrain the thermal pressure of 
the hot ionized component of the ISM at = 3.5 kpc, or in- 
deed, anywhere in the Galactic disk (see Cox 2005, for a sum- 
mary). We estimate that at the Sun's position in the Galaxy, 
Pg,Q/kB ^ 7000 cm-3 K (Cox 2005, in his §4.2). This wind 
is launched much closer to the Galactic center, so that value 
must be extrapolated to R = 3.5 kpc. Wolfire et al. (2003) 
show that the radial pressure scale-length is of order 5.5 kpc 
(notably, for the neutral gas component), so to estimate 
Pg.oiR = 3.5 kpc), we must multiply Pg ^/k^ by e'. Therefore, 
Pg.a/kYi - 1.9 X 10'* cm-3 K, ov Pg,a - 2.6 x lO''^ dyne cm'^. 

Estimates for the gas density are also difficult for R ^ 
3 .5 kpc, but for a simple estimate, if we take J ~ 10^ K for the 
hot ionized medium (that this wind would be launched from), 
and use the above pressure estimate, we find no ^ 1-8 x 10"'^, 
or p ^ 1 .9 X lO"^*" g cm"^. 

The cosmic -ray density can be estimated from the syn- 
chrotron emissivity; much of the work for this has already 
been coUected by Ferriere (2001). We used her Equations 10 
and 11, first duplicating her Figure 7 to check our own imple- 
mentation; we then use those verified equations to calculate 
both Pc and /"magnetic (^ud hence B). We find P^R = 3.5 kpc, z = 
0) = 3.1 X 10-12 jjyjje cm-2 (= 2.2 x 10^ cm'^ K). Similarly, 
/"magnetic = 2.4 X lO'^^ dyuc cm^^ (= 1 .7 X lO'i cm-3 K), which 
implies B = 7.8/iG. 

We start with a = 2 (as in BMV91), but with Zbreak = 4.5 kpc, 
as many of our early wind models preferred Zbreak of that order, 

as opposed to Zbreak = 

15inBMV91. 

2.2. Integrating Through the Critical Point 

Like many other similar wind models, this wind equation 
contains a critical point (CP) where v = c* . Such critical-point 
equations are usually solved by integrating from the singu- 
lar point to the boundaries (here at z = 1 pc and z = 1 Mpc). 
The dependence of the wind on cosmic-ray pressure and gas 
pressure mean that, in order to follow this usual procedure, 
estimates of those pressures would have to be known a pri- 
ori at the critical point position. We choose (as did BMV91) 
to estimate the pressures at the launching point of the wind, 
searching for a vq which defines the wind that threads the crit- 
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ical point. We checked our implementation by recalculating 
models presented in Table 1 of BMV91. 

This concludes our definition of the cosmic-ray and ther- 
mally driven wind model that we will compare to the soft X- 
ray data. The next section therefore returns to the X-ray data 
briefly mentioned in § 1 and reports on our search for a wind 
model to fit those observations. 

3. COMPARING THE MODEL TO ROSAT 
OBSERVATIONS 

One of the main goals of this work is to investigate whether 
a wind model can reproduce the diffuse soft X-ray back- 
ground observations most recently presented in Snowden 
et al. (1997). In this section, we describe the algorithm we 
developed to calculate that emission. 

3.1. Calculating the Average Observed X-ray Emission 

We will compare the models against the longitude-averaged 
X-ray emission as a function of latitude. We retrieved the dif- 
fuse X-ray emission maps from http://www.xraY.mpe. 
mpg . de/rosat/ survey/ sxrb/12/ f its . html at the "X- 
Ray Astronomy" page at the Max-Planck-Institutftir extrater- 
restrische Physik. These maps are the full 12' resolution 
maps; the emission maps are then averaged over the longitude 
range from -30° < / < +30°, binned 1° in latitude, and plot- 
ted vs latitude for the southern Galactic hemisphere only. The 
Northern Galactic hemisphere contains other features (such 
as the North Polar Spur) that make model comparisons there 
much less clear Note again that there is a slight asymmetry in 
the observed X-ray emission towards positive Galactic longi- 
tudes (towards the right on Fig. 1); we are not modeUng that 
asymmetry here. 

The resultant longitude averaged, observed emission is 
shown by the diamonds in Figure 5 (for the ROSAT R4 band, 
centered at approximately 0.65 keV) and Figure 6 (for the R5 
band, centered at approximately 0.85 keV). Error bars (from 
the original data files) are given by the vertical lines within 
the diamond-indicated data-points; those error bars assume 
only statistical errors in the measured X-ray flux. Systematic 
uncertainties in the foreground and background components 
(particularly the stellar contribution) may be larger. 

We restrict ourselves to the R4 and R5 bands as the emis- 
sion in those bands comes primarily from oxygen emission 
lines; higher energy bands (near 1 .5 keV, for instance) may 
depend more strongly on metallicity, as magnesium and sil- 
icon emission lines begin to dominate at those energies. In 
addition, at higher energies, the contribution from the stellar 
background becomes much more prominent and that back- 
ground is not well understood (see §3.3). 

3.2. Calculating the Wind's X-ray Emission 

The wind model gives n{z) and T{z) along ID streamlines. 
As mentioned in § 2.1, the observed large-scale emission is 
simulated with this wind model by "replicating" the wind so- 
lution both in radius (from 1.5 to 4.5 kpc) and in azimuth (see 
Fig. 4). 

First, to calculate the emissivity per emission measure, the 
emission codes "ATOMDB" and "APEC" (Smith et al. 2001) 
were used to generate spectra for a range of temperatures^ 

' As noted on their website (http://cxc.harvard.edu/ 
atomdb/), APEC is not complete below about 0.25 keV, but that will 
not greatly affect the predictions here, since the relevant bands for this work 
are at higher energies. 



from lO^ *" to 10^-^ K in steps of 10° ' K. These codes use 
the metallicities of Anders & Grevesse (1989). Each APEC- 
generated spectrum was then folded through the ROSAT re- 
sponse matrix^, with the responses of the individual channels 
summed into groups to represent the ROSAT passbands R4 
and R5 (Snowden et al. 1995)'. Applying these general mod- 
els to the wind, the resulting matrix of emissivity per emission 
measure vs temperature for each band is then interpolated to 
calculate the emissivity at each point in the wind model. We 
find that the newer APEC-derived models yield a maximum 
of 50% more emission in the M-band (the combined ROSAT 
R45 band) than the Raymond & Smith (1977)-derived models 
of Almy et al. (2000) near T = 4 x 10*' K, but the differences 
are only of order 20% near T = 2 x 10* K. 

The emission measure is calculated by simply sunmiing 
n^A/ along lines-of-sight through the wind model. This 
model emission is then corrected for absorption by applying 
(as a foreground absorption screen) the Nh data of Dickey 
& Lockman (1990) for each line of sight. The resultant wind 
emission for both the R4 and R5 bands is shown as the dashed 
line in Figures 5 and 6, respectively. 

We note briefly that we have checked that the gas in these 
winds is in equilibrium throughout the region where X-ray 
emission is important. This has been verified with the non- 
equilibrium cooling code of Benjamin et al. (2001). For in- 
stance, for the best- fit wind model presented in §3.4, we have 
found that non-equilibrium calculations yield only 1 % differ- 
ences in the population of fully-ionized oxygen vs. equilib- 
rium at z = 2 kpc, whereas for less ionized states of oxy- 
gen, non-equilibrium calculations show 1% deviations only 
beyond z ~ 4.5 kpc. Since we are interested in the X-ray emis- 
sion at z < 2 kpc, we retain simple collisional-equilibrium 
models. The previous Galactic outflow model of Breitschw- 
erdt & Schmutzler (1999), also of a hybrid cosmic-ray and 
thermal-gas pressure driven wind, relied on non-equilibrium 
effects (Breitschwerdt & Schmutzler 1994) to model the full- 
sky ROSAT emission in the all of the observed bands; the 
present model concentrates solely on the excess towards the 
Galactic Center, and for reproducing these observations and 
for this wind model, we find that coUisional equiUbrium dom- 
inates. Non-equilibrium effects may certainly be important at 
very large heights (z ^ 20 kpc) and correspondingly cool tem- 
peratures, but we do not find those effects to be significant in 
reproducing the ROSAT R4 & R5 emission. 

3.3. Calculating the Background Emission 

As first considered by Almy et al. (2000), various back- 
ground sources also contribute emission to the soft, diffuse 
X-ray observations. The stellar and extragalactic backgrounds 
listed below were used in Almy et al. (2000), and were both 
folded through both the Galactic Nh map and through the 
ROSAT response matrix. 

The extragalactic emission is calculated using 
the power-law given in Hasinger et al. (1993): 

* We used a file pspcc_gainl_25 6 . rsp, downloaded on Feb. 22, 
2006 from the "X-ray Background Tool" on the HEASARC website at 

http : //heasarc . gsf c . nasa . gov/docs /tools . html 

^ NB: We found that the channel boundaries reported in Snowden et al. 
(1995), when used to sum up the channel-by-channel responses to create the 
R4 and R5 band responses, did not recreate the response matrices plotted in 
Snowden et al. (1997) and Almy et al. (2000). The only way to reproduce 
the previous response matrices was to subtract 6 channels from the channel 
boundaries in Snowden et al. (1997). In order to compare the wind model to 
the data, it is essential that we use as similar a response matrix as possible to 
that used in Snowden et al. (1997), so we apply this channel offset 
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7.8 (E/lkeWf '^^' keV cnT^ s'^ sr^ keV"! for £ < 1 keV. 
This averaged extragalactic background is shown as the 
dot-dot-dot-dashed line in Figures 5 and 6. 

The stellar background used is the model of Schmitt & 
Snowden (1990). This background model is shown as the 
long-dashed line in Figures 5 and 6. There are significant un- 
certainties with this stellar background model; Kashyap et al. 
(1992) claim that there are errors of order a factor of three 
in Schmitt & Snowden (1990) due to the assumed luminosity 
function; unfortunately, however, Kashyap et al. (1992) do not 
model the stellar emission towards the Galactic Center. We re- 
tain the stellar model of Schmitt & Snowden (1990) because it 
addresses the stellar background in that area. The uncertainty 
in this background represents another reason for the present 
work to not address the 1.5 keV emission, which would be 
even more strongly contaminated by coronal emission; emis- 
sion in the R4 and R5 bands is not as strongly affected by the 
stellar background. 

AH of the above extragalactic and background emission, 
added together, do not fuUy account for emission at high lat- 
itudes (b < —60°). The wind model emission is also insignif- 
icant at this height. To account for this excess emission, we 
add in an isotropic component (the dotted line in Figs. 5 and 
6). This emission may be due to poorly-understood disk and 
halo emission that provides a rather uniform and apparently 
thermal-emission background at high latitudes and towards 
the Galactic anticenter This component is simply fitted (for 
each model fit to the data) to account for the emission in the 
latitude range -90° <b< -60°. 

3.3.1. Refitting the Static Polytrope Model 

The best-fit wind model will be found by comparison with 
the observations, below (§3.4), but it is also instructive to 
compare the wind model to the static polytrope model of 
Almy et al. (2000). Having updated several steps in the data 
and model analysis procedures, we re-fit the polytrope model 
to the data to ensure that both models are given fair consid- 
eration. We therefore adapted the wind-model routines and 
parameter- search codes to produce new polytrope models and 
find the best polytrope fit, again using 7 = 5/3 as in Almy et al. 
(2000). We produced a grid of 101 x 101 polytrope mod- 
els, stepping logarithmically between Pq/^b from 1.8 x 10"* 
to 1.8 X 10*^ cm-3 K and through k values of 1.45 x 10^' to 
1.45 X lO''^ cm^ g"^^^ s"'^. This parameter survey confirmed 
that the Almy et al. (2000) values for those polytrope con- 
stants remains the best-fit: we find fo/^B = 1-8 x 10^ cm"^ K 
and A: = 1.45 X lO^^ cm'' g-^/^ ^-2 j^a^ model, added to the 
background sources already considered, yields the polytrope 
models shown by the dot-dashed lines in Figures 5 and 6. 

3.4. Fitting the Wind to the Observed Emission 

The longitude-averaged wind model's X-ray emission is 
added to the various background components, and the sum 
is given as the sohd black fine in Figures 5 and 6. For each 
attempted model fit (see below), is calculated by consider- 
ing the full range of latitudes plotted, weighted by the errors 
shown in the data points. The reduced values are not close 
to unity, but we retain as a relative figure-of-merit for com- 
paring a wide range of models. 

To find the best-fit model, an initial parameter survey was 
carried out over two orders of magnitude in po and Pg o (cen- 
tered near the initial parameters given in Table 1), and a factor 
of 4 in Zbreak; the parameter ranges of po and Pgfi were grid- 



TABLE 2 
"Best-Fit" Wind Parameters 



Parameter Value Fixed? 

Pg,oAB 2.0 X 10'' K Varied 

no 6.9 X 10"^ Varied 

Zbreak 5.2 kpc Varied 

Rq Range [Galactocentric] 1.5 to 4.5 kpc Fixed 

Pc,o/*B 2.2 X 10* K cm-3 Fixed 

Bo 7.8 nG Fixed 

a 2.0 Fixed 



ded into 35 steps; while Zbreak was gridded into 5 bins. Af- 
ter this more general survey found the approximate location 
of the best-fit model, higher resolution sampling was carried 
out, leading to steps of ~ 1 .3% in po ™d Pg o, and ^ 6% steps 
in Zbreak near the \^ minimum. The best fit model parameters 
are given in Table 2 with X-ray emission shown in bands R4 
and R5 in Figures 5 and 6, respectively. The best-fit model's 
position in space is fairly well constrained, as shown by 
the green ellipse in Fig. 10. Notably, at fixed Zbreak, the fits for 
the individual ROSAT emission bands are very close together; 
the best-fit R5 model is identical to the joint best-fit R4 (& R5 
model, and the best-fit R4 model differs by only 9% in Pg^. 

Interestingly, the Pg.o value obtained is quite similar to the 
ISM value that was initially estimated. Granted, our initial, 
pre-fit estimate of Pg^ (in Table 1) could have uncertainties 
of at least a factor of two, but it is somewhat satisfying that 
the wind's required pressure is relatively close to the nominal 
ISM thermal pressure at the launch position. Meanwhile, the 
density at the base of the wind is 1/3 of our estimate of the 
ISM density. 

It is important to note also that P^-.o ^ ^g.o, and so in this 
best-fit model, cosmic rays are an important component in 
driving a wind from the Milky Way. In fact, in the best-fit 
wind model, P^. o is very slightly greater than Pg Q (but only 
by 10%). Of course, our value of PcX) is set from synchrotron 
measurements (Ferriere 2001): the diffuse X-ray emission by 
itself can place only fairly weak constraints on the cosmic-ray 
pressure, and only then because the damping of cosmic-ray 
generated waves increases the temperature of the gas slightly 
(in the best-fit model, Tmax ^ l-157b), or leads to a wind that 
cannot escape the MiUcy Way. Judging from an increase of 
X^ by ~ 50%, the X-ray data constrains Pcp to he between 
~ lO"' cm~^ K and ~ 6 x 10"' cm"-' K. (For no constant, the 
best fit Pc.o is 2.7 x 10^ cm~^ K, only 20% away from our 
assumed Pe.o-) 

The parameter Zbreak is also important to fitting the ROSAT 
data. Moving Zbreak directly impacts the location of the critical 
point, and therefore the fall-off of X-ray intensity with height. 
Having no a priori constraint for the value for Zbreak for this 
geometry, we allowed the parameter to float, finding the best 
fit value of Zbreak = 5.2 kpc. Again, judging from an increase 
in the minimum x^ by ~ 50%, Zbreak is approximately con- 
strained to lie between 4.1 kpc and 6.5 kpc. To help guide 
intuition: as Zbreak decreases, the critical point position de- 
creases, the mass outflow rate increases, and the total energy 
required increases. 

The comparison between the wind model, static polytrope 
model, and data shows that the wind is an improved fit to the 
data. The x^ values for bands R4 and R5 are smaller by a 
factor of 2.1 and 2.3 for the wind models compared to the 
static polytrope model. The x^ values are themselves very 
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Best-Fit Wind Model + Bockgrounds 

Best-Fit Static Atmosphere + Backgrounds 

O Snowden et al. (1997) Data 
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Fig. 5. — ROSAT R4-band emission (centered on ~ 0.65 keV) from the 
best fit wind model (Table 2) compared to the longitude-averaged diffuse X- 
ray emission from Snowden et al. (1997). The ROSAT data points are plotted 
as diamonds, with vertical lines representing the error bars; the eiTor bars 
are of very similar size to the plotting symbols. In the R4 band, the wind 
and static-polytrope models both fit the data reasonably well, although the 
for the wind is 2.1 times smaller than that for the static-polytrope. Still, 
systematic deviations dominate: x?, = 19.0 for the best-fit wind model in the 
R4 band. 
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Fig. 6. — ROSAT R5-band emission (centered on ~ 0.85 keV) from the 
best fit wind model (Table 2) compared to the longitude-averaged diffuse X- 
ray emission from Snowden et al. (1997). The ROSAT data points are plotted 
as diamonds, with vertical lines representing the error bars; the error bars 
are of very similar size to the plotting symbols. In the R5 band, the wind 
models fits much better than the static-polytrope model, with a difference of 
2.3 in x^- As for band R4, is not neai' unity, with a value of 48.9; again, 
systematic deviations from the model dominate. 



high (~ 5.9 X lO"' for the best fit model), showing that there 
are still many deviations between the data and wind model, 
but an examination of Figures 5 and 6 show that this simple 
wind model is at least a reasonable fit to the data, and there- 
fore merits continued detailed consideration. 

3.5. The Best-Fit Wind Model 

Before moving on to consider the range of models produced 
in our parameter surveys, it will help to consider the best- 
fit model in some detail. In Figure 7, we show the trends 
in velocity vs height for the best-fit wind model. The solid 



line represents the velocity curve for outflowing gas in the 
wind; like all thermal winds, it starts at an initial velocity 
less than the sound speed, accelerates through the critical 
point (where v = c*) and asymptotes at a velocity of order 
the escape velocity. Interestingly, in this model, the posi- 
tion of the critical point is zc? =2.4 kpc, near the edge of 
the region of observed emission. Also, this wind requires 
an initial velocity of vq = 173 km s"'; in comparison, at the 
base of the wind, c,,o = 198 km s'^ (from 7b = 2.9 x 10^ K) 
while c,_o = 251 km s"'. The terminal velocity of the wind is 
Voo ^ 760 km s~'. 
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Fig. 7. — Velocity vs. height in a fiducial wind model. The solid line rep- 
resents the wind velocity, the dashed line represents c, , the composite sound 
velocity, and the dot-dashed line shows the change in the Alfven velocity with 
height. This velocity curve shows the rather standard increase in velocity of 
a pressure-driven wind, rising from the relatively low I'o, through the critical 
point at V = c» , and accelerating on to i' ^ Voc . 



Figure 8 shows how the gas and cosmic-ray pressures com- 
pare, and how they each change with height in the wind. The 
cosmic-ray pressure and gas pressure are nearly equal at the 
base of the wind, but at large scales the cosmic -ray pressure 
drops off less quickly than thermal pressure, as we have as- 
sumed that 7c = 4/3 and 7^ = 5/3 (for cosmic rays and gas, 
respectively). The large-scale importance of cosmic rays will 
be investigated further, below. 

We present the changes in density and temperature with 
height in Figure 9. The temperature shows an increase at 
around z ^ 2-3 kpc due to the damping of the cosmic -ray 
generated waves. The density simply falls off with height 
as expected from mass conservation; the density drops off so 
quickly that the increase in temperature due to wave damping 
does not yield an increase in gas pressure at large height, al- 
though the gas pressure does drop off more slowly with height 
than it would in the absence of cosmic -ray wave damping. 

The importance of this wind for the Galaxy is largely 
summed up in the mass outflow rate: for this model, M = 
2.1 M0 yr"'. This is a very large mass outflow rate, but is in 
the same range as the inferred mass inflow rate to the Galaxy 
of r^ 1 - 10 M0 yr"' (Bregman 1999). Even with that infalHng 
gas, this wind would have important implications for metallic- 
ity gradients in the Milky Way. However, we note that while 
that mass outflow rate is well-constrained within the context 
of this simple model, there are factors which we do not con- 
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Fig. 8. — Pressure vs. height in a fiducial wind model. The solid line shows 
the thermal pressure and the dashed hne shows the cosmic-ray pressure. Be- 
cause 7r = 4/3 and 7g = 5/3 in our model, the cosmic-ray pressure drops off 
more slowly with height in the wind. In addition, as cosmic-ray driving is 
somewhat less efficient than thermal driving (see §4.3), the thermal pressure 
is most significant at the base of the wind, whereas the cosmic-ray pressure 
is most significant at large (kpc) heights. 
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Fig. 9. — Normalized density and temperature vs. height in a fiducial wind 
model. The temperature is normalizes to 10* K and the density is normahzed 
to 1 0"^^ cm"^ . The temperature briefly increases to a peak at z 3 kpc because 
of the damping of cosmic-ray derived Alfven waves. The density, however, 
drops off quickly enough with height that the gas pressure remains monotonic 
(see Fig. 8), as it must for the wind to be driven to large heights. 



sider which may significantly decrease the required mass out- 
flow rate. For instance, clumping in the wind (see §5) could 
lower the mass outflow rate required to supply the observed 
emission; the "best-fit" wind here (with constant density at the 
base) would therefore represent a maximal mass-loss wind. 

3.6. Why Does the Wind Fit Better than a Static 
Atmosphere? 

Overall, one might expect the wind shown here to yield 
emission similar to a static atmosphere, as the outflow is close 
to hydrostatic equilibrium within the critical point. This is 
true, but there is one principal reason why the wind improves 



upon the previous static atmosphere models. The wind is de- 
fined to lie between 1.5 and 4.5 kpc: the outer radial bound 
does not affect the difference of fit (since the longitude range 
of the survey data used here is limited anyway), but the inner 
bound is important, and is physically motivated: the Galac- 
tic potential there makes launching a wind unlikely, except 
under extreme conditions. (There could, of course, still be X- 
ray emitting gas within the Galactic core representing a frac- 
tion of the observed emission; the inclusion of such emission 
would lower the required mass outflow rate in this outflow.) 
In contrast to the wind model, the static polytrope must, by 
definition, fill the core of the Milky Way. This leads to ex- 
tremely high temperatures (given the deep potential) and over- 
estimates of the emission near the Galactic center in order to 
have approximately correct emission temperatures elsewhere, 
although perhaps thermal conduction could compensate for 
this (Almy et al. 2000). 

Also, we note that the drop-off in temperature with distance 
in the static polytrope model, and the necessity of the cen- 
tral very high-temperature region, means that the temperature 
elsewhere in the wind is, on average, somewhat smaller than 
required to explain the emission. This can be seen in Figure 6. 

Compared against the static-atmosphere model, then, the 
wind model is the preferred fit for the diffuse soft X-ray data. 
But how do such wind models work as the fit parameters are 
varied? These questions are addressed in the following sec- 
tion. 

4. HOW GALACTIC WINDS CHANGE WITH INITIAL 
PARAMETERS: BUILDING INTUITION 

The parameter surveys that yield a best-fit model are also 
quite useful for building intuition and understanding of mixed 
cosmic-ray and thermally-driven winds. We present a few of 
the key results below. For each survey, we have varied the 
parameters around the best-fit values for the Milky Way to 
understand how the outflow would change character near the 
best-fit parameter values. 

4. 1 . Mass Outflow Rate 

Figure 10 presents the range of mass outflow rates given 
in a variety of thermally and cosmic -ray driven wind models. 
First, it is instructive to consider the envelope of winds that 
successfully escape the Galaxy's potential vs. the unshaded 
area where winds could not be launched. Towards the bottom- 
right of the plot, at low densities and high pressures, gas is 
hot enough to escape simply by virtue of CsX) > v'esc- These 
are not outflows that our code models, and so those regions 
are not filled-in on the contour plots. In fact, if the wind were 
purely driven by thermal pressure, gas with T > Thigh in the 
plot would "evaporate" in this way. 

In the same limit of a thermal-pressure dominated wind, if 
T < Tiow, the gas would not have enough energy to escape. 
Thus, the upper limit of the shaded area in Figure 10 repre- 
sents winds that are becoming too dense (and thus too cold) to 
escape the Galactic potential. The "excess" of allowed winds 
with T < 7io„ exists because of the added cosmic -ray pressure 
gradient. The cosmic-ray pressure component acts over larger 
distances than the thermal-pressure component (as 7^ < 7g), 
and helps drive the wind where a thermal wind alone would 
fail: cosmic-ray pressure thus markedly increases the param- 
eter space where viable winds may be launched. 

Looking inside the perimeter of the shaded region in Fig- 
ure 10, we see that, as expected, as we increase the tempera- 
ture of the gas (increasing Pg.o at constant no), the mass out- 
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flow rate increases. Meanwhile, increasing no at constant f^^o 
decreases the temperature of the gas, yielding less energy to 
the gas, and so produces smaller mass outflow rates. 

In more detail, though, why does Figure 10 show curva- 
ture in the contour levels? For any given value of Pg^, there 
are two values of no where a given mass outflow rate can be 
achieved. To understand this, recall that the M contours are 
basically contours of novo. The initial velocity, vo, decreases 
as the base temperature decreases. One can think about this 
as follows: a decrease in temperature yields a decrease in the 
energy available to the gas at the base of the wind, and drives 
the mass outflow rate down; so, at some fixed Pg^, as the den- 
sity increases, vo must decrease. Since M oc novo, even as no 
increases, this decrease in vo leads to a decrease in M, and 
produces the curvature in the M contours. 

Our chosen (fixed) parameters can also affect the available 
parameter space of escaping winds. For instance, a > 2 is 
important for launching a wind; for a ~ 1 .0 the wind cannot 
pass through a critical point. Thus, the value of a can also 
strongly affect where winds are allowed. We retain a = 2 as 
inBMV91. 

4.2. Total Energy Flux in the Wind 

Of great interest for the production and impact of these 
winds is their total power (total rate of energy release in both 
kinetic energy and enthalpy). At the base of the wind, the en- 
thalpy dominates the kinetic power, so the total power simply 
scales as the base temperature in the wind (although plateau- 
ing at low gas pressure because of the added enthalpy in cos- 
mic rays). Approximately 75% (and ^ 90% for the best-fit 
model) of the power at the base of the wind is in enthalpy. 

We now check whether that power can be supplied by su- 
pernovae in the Galaxy. The total power required for the wind 
is ~ 3.7 X lO'*' ergs s~' . To calculate the normal SN power in 
the Galaxy, we use the SN rates per unit area for both Type 
I and Type 11 SN in Ferriere (2001) in her Equations 14 and 
16. Integrating over the area of the wind only, we find a SN 
rate of ^ 5400 Myr"', which equates to roughly 1 SN ev- 
ery 180 yrs. If each SN produces 10^' e ergs (where e is the 
fraction of SN power placed in cosmic rays and thermal gas), 
then the total SN power in the disk (below the wind only) is 
~ 1.7 X 10*^ e ergs s"^ This simple estimate therefore shows 
that the wind requires of order the normal SN rate in the disk, 
although it is a factor of ^ 2 too high. Of course, we are ap- 
plying a very simple model, and it is quite conceivable that, 
by including conduction (which would act as a heat source 
to add energy to the base of the wind), and by considering 
the effects of clumping within the wind, the outflow's energy 
requirements could be reduced (see §5). 

4.3. The Importance of Cosmic Rays 

Figure 11 plots the mass outflow rate in the wind, as in 
Figure 10, but as a function of cosmic -ray pressure and gas 
pressure. The soUd black contour Unes in the plot show lines 
of constant total pressure. Thus, moving counter-clockwise 
along those lines moves from gas-pressure dominated winds 
to cosmic-ray dominated winds. 

Before we consider the effect of cosmic rays in these par- 
ticular models, we start with some general considerations that 
will help us later on. First, as has been known for some 
time, momentum addition either before or after the critical 
point of a wind affects the outflow differently (e.g. Leer & 
Holzer 1980; Feldmeier et al. 1999). Momentum input be- 



fore the critical point results in an increase in the mass out- 
flow rate, whereas momentum input after the critical point in- 
creases the terminal velocity. In the context of these mod- 
els, the smaller adiabatic index of cosmic rays means that 
the cosmic-ray pressure decays more slowly (see Fig. 8), and 
therefore increases in the cosmic-ray pressure tend to increase 
the terminal velocity as well as the mass outflow rate (al- 
though cosmic rays are less efficient at "mass-loading" the 
wind, or setting the mass outflow rate, than gas pressure, as 
we wiU see shortly). On the other hand, increases in the gas 
pressure tend to increase the mass outflow rate, with relatively 
minor changes to the terminal velocity unless the gas pressure 
is dominant. This is an important point to keep in mind as we 
discuss the role of cosmic rays in these winds. 

Now consider the outUne of the shaded region in Figure 1 1 : 
this outline denotes the transition from wind models that pass 
through a critical point and winds that do not have a crit- 
ical point, or that fail to launch gas to large heights alto- 
gether. Again, at the high thermal pressure limit (and there- 
fore high- temperature limit, near Pg^ ~ 10^ cm~^ K), the gas 
is initially too hot and forms a more evaporative flow. At 
high cosmic-ray pressure towards the top of the map (P^.o ^ 
2 X 10^ cm~^ K), we see that cosmic-rays can launch winds 
even down to fairly low thermal pressures, if the cosmic-ray 
pressure dominates (more on this below). However, there is 
an upper limit in cosmic-ray pressure beyond which the mixed 
cosmic ray and thermal pressure can launch an outflow that 
does not require passage through a critical point. Beyond this 
limit, gas of all initial velocities can go to infinity, and no sin- 
gle wind solution is preferred. 

4.3.1. The Relative Efficiency of Cosmic-Ray Driving 

Looking at the basic structure of colored contour Unes in 
Figure 11, it is clear that generally, as we increase either the 
gas pressure or the cosmic-ray pressure, the mass outflow rate 
generally increases. This makes sense since increased pres- 
sure leads to increased energy in the gas which can help in- 
crease M. 

But, it is important to note that the colored contour outlines 
do not exactly follow the contours of total pressure. This is 
clearly seen at the high-pressure limit, where the outline of 
the shaded region does not follow the contours of total pres- 
sure. This is explained by the relative inefficiency of cosmic- 
ray pressure vs. gas pressure in driving gas from the Galactic 
midplane. An intuitive explanation of this comes from un- 
derstanding that (v + va)Vpc is the rate that the cosmic -ray 
pressure transfers energy to the flow (vVpe is the rate of work 
done on the gas moving at v, whereas v^Vpc is the rate of 
work done to generate Alfven waves via the streaming insta- 
bility). ForM^ <C 1, V <C V/i, and this rate is much smaller than 
downstream in the wind where v > va. As the low- velocity 
region is the region of mass-loading, it follows that M will 
drop if the proportion of cosmic-ray pressure increases for 
some given, constant total pressure. This explains why fol- 
lowing one contour of constant pressure (counter-clockwise) 
from high thermal pressure to high cosmic-ray pressure leads 
to a drop in mass outflow rate. 

In more mathematical detail, we can see the reason for this 
difference in efficiency by considering the drop in cosmic-ray 
pressure with density (and therefore drop with height in the 
wind, Eq. 4) in the sub-Alfvenic and super- Alfvenic regimes: 

Um ^ = (8) 
Ma^o dz 2p dz 
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Fig. 10. — Color contours of mass outflow rate in units of Mq yr"' , with temperature (in K) represented by solid black contour lines, as a function of base 
density, iiq, and gas pressure, Pg q. The shaded contour region shows those winds that pass through a critical point, as distinct from those regions where T is too 
high (Cs > Vcsc) and where T is too low, and the gas falls back. These regions of failed wind are approximately defined by either T > ri,i„i, or T > Ti„„: the 
temperature limits for a pure thermal wind with our parameters, shown here with the red, dashed-line contours. The excess of viable winds with T < T\a„ occurs 
where cosmic-ray pressure helps drive the wind even at low temperatures. The best fit model for both the "R4" and "R5" bands is shown as the gold cross; note 
that without cosmic-ray pressure, no such wind would be possible from the Milky Way. The area where < 2 ■ Xnim shown as the green elhpse. For winds in 
this survey, Rq, zo, Pc o. ^o> Zbreak. and a (see Table 2) were fixed. 



and lim ^ = 2^^. (9) 

M4-+00 dz p dz 

Since we assume 7^ = 5/3 and 7c =4/3, the factor of 2 in the 
the first equation above only exacerbates the difference be- 
tween cosmic-ray and thermal driving: 7^. is effectively 2/3 
compared to 7g = 5/3 in sub-Alfvenic gas. The ratio of the ef- 
fective 7 factors is 2.5; this will become important in the dis- 
cussion below. So, in the sub-Alfvenic regime, dPc/ dz drops 
more slowly than dPg/dz, which means that less momentum is 
imparted to the gas; hence, we may conclude that cosmic rays 
drive gas less efficiently when the wind is sub-Alfvenic. If the 
winds were everywhere (in Fig. 11) launched sub-Alfvenic, 
then the cosmic-ray driven winds would require pressures 
about a factor of 2.5 higher than the thermally-driven winds. 

This complication is that, at high P^.o ~ 10^ cm"^ K, the 
winds are somewhat super-Alfvenic at the base, whereas at 
lower Pc.Q, the winds are Alfvenic to sub-Alfvenic (Fig. 1 1 is 
for fixed no, so the colored contours on this plot are curves of 
vo, which therefore decreases as Pc.o decreases). Winds that 
are launched with super-Alfvenic velocities will have more 
efficient cosmic-ray driving. So, some of the difference in 
required launching pressures between cosmic -ray and ther- 
mally dominated winds are also due to a transition from sub- 
Alfvenic to somewhat super-Alfvenic which compensates for 
the fact that the winds are not strictly launched in the Ma ^ 1 
ox Ma ^ 1 regimes. 

In understanding the relative efficiency of cosmic-ray driv- 
ing, we have only explained the basic trends in M with to- 
tal pressure and the high-pressure limit of Figure 1 1 . We 



now turn to the low-pressure Umit, where winds also cannot 
emerge from the Galactic midplane. 

4.3.2. The Impact ofAlfven Wave Damping 

The relative efficiencies of cosmic-ray and thermal pres- 
sure driving do not explain the lack of winds at low thermal 
pressure (P^ < lO'* cm"^ K) and low to intermediate cosmic- 
ray pressure (Pc,0 < 10^ cm"-* K). At the very lowest total 
pressures (Ptotai ^10"* cm"-' K), winds cannot be launched be- 
cause the pressure is simply too low (cosmic -ray driven winds 
would be possible for lower initial densities, however). But at 
higher cosmic -ray pressures, winds cannot be launched be- 
cause Alfven wave damping and the resultant heating dom- 
inates. For intermediate P^.o and low Pgo, the damping of 
Alfven waves can actually lead to a large increase in tempera- 
ture and therefore an increase in Pg with height; this pressure 
increase can cause the wind to stall (BMV91). So, in this 
regime, the Alfven wave damping actually "pressurizes" the 
Galactic disk and prevents a wind from being launched from 
the midplane. We have verified this result by running a pa- 
rameter survey identical to that shown in Figure 1 1 , but with 
wave dissipation turned off; as expected, the lack of dissipa- 
tion allows winds to form with P^ o < lO"* cm"^ K and with 
Pc,0 < 10^ cm"^ K. For winds with higher Pc.o and low Pg o, 
the cosmic-ray pressure dominates to such an extent that the 
wave damping does not hamper wind driving. (Of course, all 
of these constraints on pressure components are only strictly 
valid for launching from the Milky Way's midplane, as we 
have assumed.) 
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4.3.3. Launching Winds at Low Pc, a 

We now understand most of Figure 1 1 except for the small 
gap in wind models at low cosmic -ray pressure and interme- 
diate gas pressure (10"* cm"-^ K < Pgf^ ^ 4 x 10"' cm""* K, 
Pcfi ^5x10^ cm"^ K). In this region, the winds are not so 
strongly dominated by either cosmic-ray pressure or gas pres- 
sure, Alfven wave damping will not be important (because 
of the relatively low Pc), and the two must operate in con- 
cert. The unshaded region defined above is essentially the 
region where the inequality in the pressures leads to diffi- 
culty in self-consistently launching a wind (which gas pres- 
sure is especially good at; see above) and driving it to infinity 
(where cosmic-ray pressure starts to dominate). If the ther- 
mal wind attempts to load too much mass at the base of the 
wind (more than the very low cosmic-ray pressure can han- 
dle), the wind fails at large radii. If the thermal wind loads the 
wind with a fairly low mass outflow rate, the very low cosmic- 
ray pressure can continue to loft the gas to large distances 
after the gas pressure decays away. This is corroborated by 
the very low mass outflow rate associated with the "promon- 
tory" of low mass outflow rate near Pg^ = 9 x 10^ cm"^ K and 
Pcfi = 2 X 10^ cm-3 K. 

5. CONCLUSIONS 

We have implemented a simplified cosmic ray- and 
thermally-driven wind and have used it to try to explain the 
soft, diffuse X-ray emission seen towards the Galactic Center. 
We find that such a wind can indeed match the observed aver- 
aged X-ray emission quite well, and in fact fits demonstrably 
better than the static polytrope model of Almy et al. (2000). 
It is important to note that this wind is approximately equally 
powered by both cosmic rays and thermal pressure: cosmic 
rays are important in helping this relatively cool wind escape 
from the Galactic potential. It is also quite interesting that this 
wind does not require excessive thermal or cosmic -ray pres- 
sures (both pressures are not extreme compared to what has 
been estimated for the inner Milky Way), nor does this simple 
model require much more energy than the standard inferred 
supernova rate implies. Taking this result at face value, such 
a wind would be very important to the "ecology" of the Milky 
Way due to the high mass loss rate of 2 Mq yr"'. In addi- 
tion, such a wind would also play an important role in remov- 
ing angular momentum from matter in the Galactic disk and 
allowing matter to move radiaUy inward (Zirakashvili et al. 
1996). At the least, this shows that such wind models should 
be considered further for the Milky Way; they may be able 
to explain at least a substantial fraction of the observed soft 
X-ray emission. 

Further, as other researchers have already shown (Bre- 
itschwerdt et al. 2002), such winds can also be used to explain 
the unexpectedly slow rise in 7-ray emission towards the cen- 
ter of the Galaxy; this work therefore gives independent sup- 
port to the Galactic wind hypothesized in Breitschwerdt et al. 
(2002). We have not yet calculated the effect of the best- fit 
wind model on the cosmic-ray distribution, but a simple ap- 
proximation will investigate if this wind is removing cosmic 
rays at too high a rate. We reason as follows. The best-fit 
wind model shown here has a cosmic-ray advective timescale 
of ~ 4.3 X 10^ years. Therefore, supernova in the disk must 
resupply the cosmic-ray pressure on that timescale. We know 
that the approximate total supernova energy in the disk (from 
observations, see §4.2) is ~ 1.7 x 10"*' ergs s"'. So, if some 
fraction of this energy, ccr, is given to cosmic rays, and dis- 



tributed over the volume occupied by hot gas where the wind 
is launched (over the Galactocentric radius range of 1.5 to 
4.5 kpc, and a height range of ±2 kpc, to be conservative), that 
energy density should be similar to the cosmic-ray pressure 
required to launch the wind. Calculating the resultant buildup 
of Pc over 4.3 x 10^ years at the observed SN rate, we find 
Pc - ecR • 6.9 x 10-'^ dyne cm'^ or - ecR • 5.0 x lO" K cm"l 
This is actually of order the P^ that the best-fit model requires 
(see Table 2), although it would require ecR ^ 0.6 to dupli- 
cate the best-fit Pc, which is relatively high. Still, this simple, 
rather conservative calculation shows that the high M wind 
shown here does not remove cosmic rays much more quickly 
than they can be replenished by the normal SN rate in the 
Galaxy, although the removal rate of cosmic rays is certainly 
non-negligible, and would affect the density of cosmic rays 
towards the Galactic center. Of course, a more detailed calcu- 
lation is required (with a more detailed wind model), but this 
again shows that the best-fit wind model is at least feasible, 
and would have a significant but not destructive effect on the 
Galaxy's cosmic-ray density. 

5.1. Future Improvements 

More detailed models are clearly needed; there are a few 
concerns about the current model that could be addressed with 
more reaUstic wind models. For instance, we have assumed 
uniform density at the base of the wind over the area of the 
disk from R=\.5 kpc to 4.5 kpc, which leads to a mass out- 
flow rate of order 2 Mq yr"' . This seems quite high, but in 
the context of a more detailed model with variations in den- 
sity within the wind, we might expect that the weighting 
of emissivity would favor overdensities, and allow an inho- 
mogeneous wind to better reproduce the observations with a 
smaller mass outflow rate. In addition, it is possible that other 
effects limit the gas to velocities below those in this simple 
model; drag effects may slow down the wind (e.g., Everett & 
Murray 2007), and could lead to some of the gas forming part 
of the Galactic fountain (Bregman 1999). On the other hand, 
turbulence may be an important additional source of energy 
for the wind, but we have not included such an input in this 
work. Also, the effects of distributed mass loading, which 
could be relatively easily incorporated into this model, have 
been ignored so far; such mass loading would be very impor- 
tant to the emission properties of the wind, and the potential 
observability of such winds in many galaxies. Finally, note 
that we have used Anders & Grevesse (1989) abundances; if 
the wind starts out with super-solar abundances (particularly 
in oxygen), a smaller mass outflow rate would be required. 

The cosmic ray physics in this wind is still quite simple. 
As mentioned previously, we have assumed zero diffusivity 
of the cosmic rays throughout the wind (c.f., Breitschwerdt 
et al. 1993). We have also neglected thermal conductivity 
(our calculations show that conductivity may be important for 
z < 350 pc). Including both of these effects would be essen- 
tial to future progress for this wind model: as mentioned pre- 
viously (see §4.2), including the effects of conduction may 
help lower the energy requirements in the wind, bringing the 
total power of the wind closer to that of the total inferred 
supernovae power in the Galactic disk below the wind. We 
also note that we are using a model of the gravitational po- 
tential that we know overestimates the potential in the Milky 
Way (see §2). Adopting the more detailed potential model 
of Dehnen & Binney (1998) may allow a lower total-energy 
wind to duplicate the soft X-ray observations. Finally, we 
have also assumed the hydrodynamic model of McKenzie & 
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Fig. 1 1. — Mass loss in the wind vs. clianges in the initial cosmic-ray pressure and gas pressure. The black contour lines represent lines of constant total 
pressure. The winds in the shaded region of the contour plot represent mass loss in outflows that pass through a critical point; the unshaded regions show locations 
of parameter space where winds cannot be launched (where the total pressure is too low) or "evaporative" flows that do not pass through a critical point (at very 
high total pressure). For winds in this survey, Rq, zo, «o. Bo, ^■break' ^ (sec Table 2) were fixed. 



Webb (1984) and Breitschwerdt et al. (1991) for the interac- 
tion of cosmic rays with Alfven waves and the gas. To fur- 
ther examine this model, we will next consider the effects of 
higher cosmic-ray fluxes (Zweibel 2003) and apply it in other 
settings. 

5.2. Future Tests 

How can we further test this model? In analogy to early 
studies of the solar wind, this outflow may impact clouds in 
the vicinity of the galaxy, perhaps causing "comet-tail" ex- 
tensions to high velocity clouds above the plane of the Milky 
Way. The formation of such "tails" would depend on the ve- 
locity of the wind. This has been studied in some detail be- 
fore (Benjamin & Cox 2002), but should be reconsidered in 
the context of the predictions of these winds. In addition, it 
may be possible to study the kinematic impact of this wind 
on, for instance, the Magellanic Stream (A. Burkert, personal 
communication). 

Another way to test the model would be to compare the 
wind with absorption columns and emission spectra towards 
the center of the Galaxy. Concerning absorption measure- 
ments, recent Chandra observations (Futamoto et al. 2004) 
towards the low-mass X-ray binary 4U 1820-303 show sig- 
nificant columns in O Vll , O Vlll , and Ne IX . As this X-ray 
binary is located within ten degrees of the Galactic center, at 
a distance of approximately 7.6 kpc, it seems an ideal target. 
However, similar absorption columns are found towards ob- 
jects on sightlines that do not intercept the base of the wind; 
for instance, Mrk 421 shows similar absorption columns (Fu- 
tamoto et al. 2004), but is far out of the plane. This leads us 
to conclude that the observed absorption is somewhat local 
to the Sun's position, and, as such, this absorption does not 



constrain the wind model. However, recent Suzaku measure- 
ments of emission at various latitudes along different sight- 
lines towards the Galactic center (Rocks, 2008, in prepara- 
tion) may help constrain the properties of the wind. 

Finally, these winds, including their important cosmic -ray 
component, will also emit synchrotron radiation. We are 
now calculating the synchrotron emission expected from these 
wind models (Schiller et al., 2008, in preparation). This will 
allow exploration of the wind's synchrotron emission as com- 
pared to recent models of Galactic synchrotron which begin to 
map the three-dimensional cosmic-ray emissivity (Nord et al. 
2006). 

Looking at a wider field of application, such wind models 
may also be quite important in application to starburst galax- 
ies (e.g., Gaflagher & Smith 2005; Socrates et al. 2006) and 
dwarf galaxies. The general applicability of these kinds of 
models to starbursting galaxies has been shown by Breitschw- 
erdt (2003) in fitting cosmic-ray and thermally driven wind 
models to NGC 3079. 
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